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Abstract 

Recent improvements in the radiation transport code 
HZETRN/BRYNTRN and galactic cosmic ray environmental 
model have provided an opportunity to investigate the effects 
of target fragmentation on estimates of single event upset 
(SEU) rates for spacecraft memory devices. Since target 
fragments are mostly of very low energy, an SEU prediction 
model has been derived in terms of particle energy rather than 
linear energy transfer (LET) to account for nonlinear 
relationship between range and energy. Predictions are made 
for SEU rates observed on two Shuttle flights, each at low and 
high inclination orbit. Corrections due to track structure 
effects are made for both high energy ions with track structure 
larger than device sensitive volume and for low energy ions 
with dense track where charge recombination is important. 
Results indicate contributions from target fragments are 
relatively important at large shield depths (or any thick 
structure material) and at low inclination orbit. Consequently, 
a more consistent set of predictions for upset rates observed in 
these two flights is reached when compared to an earlier 
analysis with CREME model. It is also observed that the 
errors produced by assuming linear relationship in range and 
energy in the earlier analysis have fortuitously canceled out 
the errors for not considering target fragmentation and track 
structure effects. 

I. INTRODUCTION 

The demands for reducing the size and power requirements 
of future spacecraft leads to designs which efficiently process 
large quantities of data onboard to reduce the amount of 
telemetry. Such an approach leads to designs utilizing high 
density digital devices which are more susceptible to single 
event upsets (SEU). Accurate prediction of SEU rates for 
spacecraft computers in orbital environments is needed for 


cost-effective system designs. The accuracy is determined by 
several factors, such as device interaction models, accelerator 
measurements and interpretation of measured data, 
environmental models external to the spacecraft, and transport 
calculations for modified environments internal to the shield. 
Although progress has occurred in these areas, as described 
and summarized by Petersen et al. [1], uncertainties in existing 
external and internal environmental models are still very large. 
Recently, improvements have been made by Badhwar and 
O'Neill in predicting the free-space galactic cosmic ray (GCR) 
energy spectrum [2]. Their model was tested against the LET 
spectra measured on several Space Shuttle flights [3], A 
detailed comparison of the measurement with tissue 
equivalent proportional counters (TEPC) with predictions 
obtained using the model [2] and the radiation transport code 
HZETRN/BRYNTRN [4-7] shows overall agreement to 
within ±15 percent for the absorbed and equivalent dose but 
with some discrepancy in LET spectra. The underprediction 
at lower LET suggested the need for including the secondary 
pion and kaon fluxes in the code and developing an improved 
understanding of the detector response due to delta rays and 
wall effect [3]. The underprediction in the higher LET region 
requires further refinement of the energy spectra of target and 
knock-on particles. Recently, modifications have been made 
in the code to include explicit description of the energy 
spectra of heavy ion target fragments and refinements in those 
of the light ion secondaries [7,8] thus improving the 
calculation for high LET spectra, critical to SEU evaluations. 

Target fragments are mostly produced by energetic 
protons, neutrons, and some light ions in the LEO 
environment. In the primary GCR, ^He is the most abundant 
nucleus after 1 H. Since the GCR spectrum is broad with a 
peak at several hundred MeV/A with an appreciable number 
of particles out to about 10-50 GeV/A, fragmentation of the 
multiple charged ions is dominant in GCR transport. The light 
mass particles have longer ranges and great multiplicity in 



elastic nuclear events, leading to a buildup in their number in 
shielding. Estimation of upset rates due to these particles 
becomes important when GCR and trapped protons pass 
through a bulk material. Rather than relying on costly 
laboratory SEU cross section measurements for these 
particles, several investigators [9-11] have attempted in the 
past to predict proton or neutron SEU rate from low energy 
heavy-ion test data. The methods used to obtain nuclear 
reaction products were mostly time-consuming Monte Carlo 
simulations, unlike HZETRN/BRYNTRN code which 
calculates nuclear reaction products very efficiently. 

O'Neill and Badhwar recently made a preliminary analysis 
for the measured SEU rate of SRAM memories situated at the 
Space Shuttle general -purpose computers (GPC's) for several 
shuttle flights [12]. The transport calculations [12] were made 
using the CREME (Cosmic Ray Effects on Microelectronics) 
code [13], with the model from reference [2] serving as the 
input GCR spectrum. The results indicate a need for a more 
exact transport calculation which would account for the target 
fragments from the breakup of shield nuclei, particularly for 
low inclination orbits. It is anticipated that inclusion of target 
fragmentation will reduce the inconsistency in their 
comparison between low and high inclination orbits. In this 
work, we will revisit the problem for incident GCR by 
developing a SEU prediction model suitable for low energy 
target fragments and describing the necessary corrections for 
track structure effects on the available charge for upset. The 
transport calculation has been made using the newly modified 
HZETRN/BRYNTRN code which calculates very efficiently 
both projectile and target fragment energy spectra for 59 
isotopes from neutrons up to iron ions. 

II. Galactic Ion Transport Model 

The propagation of the GCR and their secondaries through 
bulk matter is described by the Boltzmann equation which in 
the straight ahead approximation is of the form [5,14] 
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where Vy denotes the range scaling parameter which is equal to 
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Zj / Aj ; where Aj and Zj are the mass and charge numbers of 

ion j, respectively. In equation (1), E represents energy 
(MeV/A), S(E) is the proton stopping power, <5j(E) is the total 
cross section, (fy (x,E) is the differential flux spectrum, and 
fjk(E,E') is the differential energy cross section for 
redistribution of particle type and energy through elastic 
scattering or nuclear reaction such as fragmentation. The 
numerical solution to equation (1) has been developed by 


Wilson et al. [14] using the method of characteristics with the 
production terms separated into projectile fragmentation and 
target fragmentation terms. 

The HZETRN code (for which BRYNTRN is a subset for 
the light particle transport) uses energy dependent nuclear 
interaction cross sections and assumes realistic energy spectra 
for low mass particles (A < 4) [7,15]. Inclusion of the energy 
dependent nuclear cross sections was shown to have a large 
effect on the highest LET components [16] and is essential to 
predicting SEU rates from GCR. For heavy ion projectile 
fragmentation, secondaries are assumed to be produced at the 
velocity of the primaries. Further details on the transport 
methods and data base are found in [5,15,16]. 

The input free-space GCR spectrum is calculated using the 
model of Badhwar and O'Neill [2]. The value of the 
deceleration potential required as input to this model was 
derived from the measured Climax neutron monitor count rate 
which varies over the 22-year solar cycle. The standard GCR 
model commonly used to predict SEU rates is the Cosmic Ray 
Effects on Microelectronics (CREME) code [13] which 
accounts only for 11 -year solar cycle modulation (represented 
by a simple sine variation) and does not distinguish between 
odd-even cycles with 22 year period. For low earth orbit, the 
geomagnetic cutoff given by Shea and Smart [17] later coded 
by Adams et al. [18] is used. The HZETRN and the model of 
Badhwar et al. have been tested with several Shuttle flight 
experiments using passive detectors (TLD and nuclear track 
detectors) as well as active detectors (tissue equivalent 
proportional counters and a charged particle spectrometer 
telescope) and appeared to be in reasonable agreement 
between model predictions and measurement overall [8] but 
some differences remain unresolved as described in Section 1 . 

We further compare the results of the HZETRN code using 
environmental model of Badhwar and O’Neill with the 
measurements of Wiegel et al. [19] using CR39 track 
detectors. Their measurements include both the contributions 
of the HZE particles (primaries and their fragments) and the 
low energy target fragments. The short range target fragments 
were identified as adjacent etch cones across the interface of 
two track foils and no entrance track at the top of the stack. 
The etching process removes 40 pm at the foil interface so 
that only those tracks of 40 pm and greater length are counted. 
In addition, many of the isotropically produced target 
fragments not counted must have greater range if they are 
oblique to the interface. We account for the etching loss of 
short range tracks resulting in reduced contributions from low 
energy fragments. Figure 1 shows the comparison of 
measured LET spectra for D1 mission (first German Spacelab 
mission) obtained by Wiegel et al. [19] with the current 
environmental and transport model. There is a vast 
improvement of HZETRN with its current environmental 
model to predict the CR39 result over the results given in 
Figure 1 of reference 19. Unlike earlier calculational models 



[19], the trends in the measurement are well represented by 
the calculations at a single depth of 20-30 g/cm-. Not least of 
these improvements is the fragmentation model improvements 
in the lower LET region. 



along various rays with a distribution of chords of length s 
given as f(s). The number of particles with chord length 
between .v and s + ds is 

dNi = n A f(s) ds ()/ ( E ) dE (4) 

The energy loss by the particle along the chord is 

e = E-R; ] [R i (E)-s] ( 5 ) 

where R t (E) is range of particle i of energy E, and R [ 1 (x) is 
the residual energy for a particle i of residual range x. (Note, 
R and R 1 are operators such that Rj 1 [/?,•(£)] = E). 
Associated with the critical energy is a critical chord s c given 
by 

*c = E-Ri 1 [W -*c] (6) 

which may be solved as (for E > 8 r ) 


Fig. 1 Comparison of measured LET spectra for D1 mission 
with HZETRN calculation with solid lines representing 
spectra behind various thicknesses of Al. 


III. Single Event Upset Model 

Direct Ionization 

The device upset results from the release of energy by ions 
to generate electron-hole pairs resulting in a current of 
sufficient magnitude to switch state in the device. It is usually 
assumed that the single parameter critical charge is 
independent as to where the free charges are formed within 
the sensitive volume; whereas in reality there is a charge 
collection efficiency and charge recombination effect (for a 
dense track) depending on the initial charge distribution. The 
charge generated is related to the energy deposit as 


s c = R 'i (E) - Ri (E-£ c ) (7) 

At the lower energy end of a typical GCR spectrum, all 
particles of energy E will cause an upset when traversing the 
sensitive volume along chords of length greater than s c . At 
the high energy end, only those particles depositing energy 
greater than e c after traversing a distance limited by the 
longest chord s max of the sensitive volume will cause an 
upset. The total number of upsets is 

dN = r max f(s)ds n A(L (E) dE (8) 
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Note that s c depends on both particle type and energy. 
The total upsets caused by all particles and energies is 


i 


n s 

max f(s) dsnMpAE) dE 
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Q (in pC) = e (in MeV) / 22.5 (2) 

It is customary to assume that an upset occurs if Q > Q c 
corresponding to a critical energy e c = 22.5 Q c . In this case, 
the upset event for a given particle traversal is reduced to 
estimating the energy loss in crossing the sensitive volume. 

In a uniform isotropic radiation environment as often 
found in a spacecraft, the number of particles entering the 
device sensitive volume is 


which may be calculated as 

poo 

A = E F(s c ) nA^i(E) dE (10) 
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where F(s c ) is the integral pathlength (chord) distribution. 

Another form of equation (10) can be obtained by first 
approximating the critical chord by a Taylor series as 


dNj= kA §i(E) dE 
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s c = £ c /L(E) + 0(£ l) (11) 


where ^ (E) dE is the omnidirectional fluence of particles of where L is the linear energy transfer ( LET ). Retention of the 
type i between energies E and E + dE and A is the total surface fj rs t term is equivalent to assuming L being constant along the 
area of the sensitive volume. The particles arrive randomly chord One may rewrite N as 
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where <J>; (L) is the integral LET spectrum. A change of 
variables as L = e r /,v results in 

N='LnAe i [ 0 °/(e c I L) <t> ; (L)dL/ L 2 (13) 
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consistent with the results used in the CREME code. The 
CREME results require validity of the Taylor series expansion 
which is only satisfied at high energies for Rf (E) » s c . 

Since most test data do not produce a sharp step function 
for the SEU cross section measurement, it is customary to 
understand this behavior in terms of distributed device 
elements with varying upset sensitivity. We follow the usual 
custom and perform an expected value calculation [20] for the 
upset rate N by integrally weighting N(e c ) with the 
normalized experimental cross section curve 

poo 

A = J e N(e c ) g(e c ) de c (14) 


emanating from the core of particle track varies widely 
according to the energy and charge of passing ions as seen in 
nuclear emulsion. At very large energy per unit mass, the 
track width is usually many microns causing energy deposit 
beyond the sensitive site (carried by secondary electrons). In 
this work, the fraction of energy deposited within the sensitive 
volume over the energy loss as given by Xapsos [21] will be 
used to correct for large particle track width. For low energy 
high-LET target fragments, the charges created along the core 
are so dense that charge recombination reduces the amount of 
available charge for upset. A simple charge generation and 
recombination model developed for the high-LET fragments is 
described below. 

The energy release by the passing ions is through 
collisional excitation of orbital electrons of atoms within the 
device to high momentum states. These electrons in turn 
diffuse from the collision site transferring significant amounts 
of the kinetic energy received from the ion to other electrons 
as the initial energy released is degraded and distributed about 
the ion path. The electron slowing down spectra are 
remarkably independent of the initiating event as the 
collisional transfer among electrons dominates the diffusion 
process. Ultimately, the energy resides in subexcitation 
electrons which thermalize and recombine with the holes 
produced in the media or are passed through the external 
circuit of the device. Recombination is dominated by Auger 
transitions and follow the kinetic equations 


where £ m i n is the minimum threshold energy and g (e c ) is the 
probability density distribution function related to the 
cumulative distribution function as represented by the 
normalized experimental cross section curve. We further 
approximate the integral over e c with summation by dividing 
the measured cross section curve into several steps thus 
yielding small surface areas. Substituting equation (13) into 
equation (14) one obtains 

N = nL Aj £ J £max F(s c ) <]>,(£) dE (15) 

j i Je v 

where Aj is the total surface area of step j and £ c j = L t ] 7 j • t 
with t being the charge collection length and L t } r j corresponds 
to the LET threshold for the /-th step. Note that if one uses 
equation (10) rather than equation (13) for equation (14), one 
obtains an SEU rate expression similar to equation (2) in 
reference [12] except that the charge collection length t does 
not change with the choice of steps but rather is given by the 
real physical size. Since equation (10) is more accurate than 
equation (13) in evaluating the low energy target fragments, 
equation (15) is used in this work. 

Track Structure Effect 
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dt 


(16) 


where we assume the electrons and hole densities are similarly 
distributed. The time development of the electron density is 

n(t) = (17) 

1 + /ZqEt 


We follow Martin et al. [22] and consider the degree 
of recombination in the device charge recombination time x 
before charge collection occurs where (Ex) - ^ ~ 5x 1 0^* /cm^ 
for 10 ps. The initial density can be found from the work of 
Cucinotta et al. [23] as 


n 0 ( r ) = D ( r ) !W 


(18) 


where 


D(r) 


2.33x0.085Z 2 10 7 r 
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1 
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The response of materials to passing ions is related to the 
amount of local energy deposited and the manner in which 
that energy is deposited. The scale of delta-ray track 


is the absorbed energy per unit volume, r max is the maximum 
electron range resulting from the ion collision, and 



W is 3.6 eV per ion pair formed. Thus, the collectable charge 
distribution (after recombination) is approximately 


n(%, r) 


np(r) 

l + n Q (r)lcx 


( 20 ) 


from which the total charge produced by the ion passage can 
be calculated. The initial charge produced per unit pathlength 
of the ion is 


poo 

Qo =J dr 


( 21 ) 


thicknesses (3, 5, 8, 12, 20, 30, 40, 55, 70, 85, and 100 g/cm 2 ) 
that cover the range of thickness for the computer location in 
the orbiter. The aluminum equivalent thickness of the orbiter 
body was determined for 968 rays uniformly distributed in all 
directions emanating from the computer location. The 
fraction for the number of rays (or fraction of solid angle) at 
these eleven thicknesses are given in Table 1. The energy 
spectra resulting from the transport calculation at each of the 
eleven depths are then used in Equation (15) to obtain SEU 
rates which in turn are weighted with the distribution given in 
Table 1 to give the combined SEU rate predicted for the 
orbiter. 


The fractional charge remaining after recombination for 
several ions is given by 

—— = \ Q 2 nrn(x,r) drj j^ 2 nrn 0 (r) dr ( 22 ) 

Figure 2 shows the results for a charge recombination time 
of 50 ps which is arbitrarily taken for illustrative purpose. 
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Fig. 2 Fractional charge remaining in ion track core as 
function of (3 (ion velocity divided by speed of light) 
for various ions of charge Z after 50 ps of 
Auger recombination. 


Table 1. Shuttle Shield Distribution 


Equivalent aluminum thickness, g/cm 2 

Fraction of solid angle 

3 

.100 

5 

.142 

8 

.086 

12 

.054 

20 

.082 

30 

.109 

40 

.113 

55 

.170 

70 

.055 

85 

.045 

100 

.044 


The experimental SEU cross section curve for 
IMS1601EPI used in previous analyses [12] is shown in 
Figure 3 and is used in equation (15) to predict error rates. 
Here we still assume the geometry of sensitive volumes to be 
rectangular parallelpiped with an assumed thickness to width 
ratio of about 0.1, yet at each step for carrying out the 
summation in equation (15) the critical energy is related to the 
LET threshold at each step with an thickness of about 1 pm. 
The degree of charge recombination assumed here is 
(Uc)“l = 5xl0 2 ®/cm 2 . 


IV. Shuttle Data Analyses 

The Space Shuttle general purpose computers were 
upgraded with IMS1601EPI SRAM to replace original 
magnetic core memories in 1991. The contributions to SEU 
rates recorded in many shuttle flights were previously 
analyzed and separated [12] into two sources: GCR and South 
Atlantic Anomaly. In this study, only the contribution from 
GCR on two shuttle flights are considered: STS-56 
in 57.1 degree inclination orbit and STS-51 in 28.5 degree. 
The input GCR spectra with geomagnetic cutoff are 
propagated through 100 g/cm 2 of aluminum shield using 
HZETRN/BRYNTRN code. Shuttle orbiter shielding is 
accounted for by obtaining energy spectra at eleven different 



Fig. 3 Measured heavy ion single event upset integral 
cross section versus effective LET for IMS1601EPI 
as given in reference 12. 




v. results and Discussion 

A calculated LET spectra for projectiles and their 
fragments, excluding target fragments other than nucleons, is 
shown in Ligure 4 for GCR integral fluxes propagated through 
aluminum shield thicknesses of 55 g/cm^ at STS-56 flight 
conditions. Also shown are the spectra with light ion 
fragments (Z < 3) added, and with all the fragments (all Z) 
added. The increase due to heavier target fragments is above 
1050 MeV/g/cm^ while increase caused by light ions is 
slightly below that. The lowest threshold for SEU seen in 
Figure 3 is 2.5 MeV/mg/cm^ (2500 MeV/g/cm^). With the 
sensitive volume thickness to width ratio of 0.1, an increase in 
the integral flux above 250 MeV/g/cm^ will affect the 
calculated SEU rate. Although the increase in flux due to 
added target fragment contributions is about one to two orders 
of magnitude larger than without in the LET region of interest, 
their contributions are greatly reduced by the track structure 
effect shown in Figure 2. 



Fig. 4 Calculated integral LET spectra at 55 g/cm- depth in 
aluminum shield exposed to GCR at STS-56 flight conditions. 

Contributions from all the target fragments added, light ion 
target fragments, and no target fragments added are 
shown separately. 

The production and accumulation of light ions as GCR 
propagate through a shield are illustrated in Figure 5 for 
STS-56 environment. The flux for primary ^Hc particles with 
a peak near 1000 MeV/A is seen to decrease with the depth in 
the shield while the secondaries from target fragmentation 
accumulate at low energies which correspond to high LET. 


The isotope ’He is shown to accumulate with low energy 
fragments. A similar situation applies to deuteron and triton 
which has less contributions to SEU rates because of lower 
stopping power than helium. 

The modified GCR fluxes due to atomic and nuclear 
reactions with the Shuttle shield material as described above 
are then used in equation (15) to predict SEU rate at a given 
shield depth. Some of the results are shown in Tables II and 
III for STS-56 (57.1 degree inclination) and STS-51 (28.5 
degree) flights, respectively. The contributions from various 
radiation components are also listed. In general, at thinner 
shield thicknesses, the rate is dominated by the projectiles and 
their fragments, but for thick shields light and heavy target 
fragments dominate. For the lower inclination orbit (STS-51), 
the ratio of SEU rates from target fragment contribution to the 
total is generally higher since there is a higher proportion of 
high energy incident GCR particles which produce target 
fragments as compared to high inclination orbit. Note that the 
representation by equivalent aluminum thickness for all the 
orbiter walls and other structure materials in current and prior 
[12] analysis is only an approximation as far as nuclear 
reaction is concerned (see, for instance, reference [27] about 
the effect from choice of material). 

The individual upset rates calculated above for each 
thickness are then weighted over the orbiter shield distribution 
shown in Table I to obtain the final SEU rates as listed in 
Table IV. Also listed are the flight data and previous works 
[12] for comparison. The observed upset rates on Shuttle 
flights have been reduced in reference [12] into separate 
contributions from GCR and SAA (South Atlantic Anomaly) 
with the GCR rates listed in Table IV. In our present analysis 
the ratio of SEU rates for two orbits is 3.8:1 versus 2.8:1 from 
flight data while the prior analysis gives 4.7:1, indicating that 
the current analysis is more consistent with the flight data 
when considering the differences in the two incident GCR 
spectra. Note that a good agreement between flight data and 
prior analysis [12] in absolute values is fortuitous since the 
analysis was based on linear relationship between range and 
energy used in their model and the charge collection length 
embedded in their model can vary indefinitely with choice of 
summation steps as discussed in Section III. Furthermore, it 
was shown in Figure 9(a) of reference [12] in comparing their 
predicted upset rates with flight data that no consistency 
between sets of low and high inclination orbits can be reached 
when they vary their device sensitive volume thickness to 
width ratio. The listed results in Table IV for their analysis is 
for their best compromised value of device sensitive volume 
thickness to width ratio. 
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Fig. 5(a) 4 He 
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Fig. 5(b) 3 He 
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Fig. 5 Production and accumulation of light ions 
at GCR propagates through aluminum shield at 
STS-56 flight conditions. 

In current analysis, there are several uncertainties 
associated with the calculated results. Obviously, the degree 
of charge recombination and charge collection length are not 
exactly known. The crude selection of these parameters is 
permissible when considering other assumptions existing in 
conventional SEU model, such as rectangular parallelpiped 
shape for sensitive volume [28]. Errors are also often 
introduced in reducing measured cross section data into a 
single curve [20], Note that the cross-section data shown in 
Figure 3 do not allow further refinement to the existing curve 
for the present analysis since the raw data were not available. 


Table II. SEU Contribution From Various Radiation 
Components Behind Several Shield Thicknesses for STS-56 
(upsets per day per computer) 


Aluminum thickness, 
g/cm 2 

0 

3 

5 

20 

55 

100 

Projectiles and 
projectile fragments 

7.08 

7.48 

6.95 

3.61 

.78 

.12 

Light target 
fragments (Z=l,2) 

0 

.70 

.73 

.91 

1.17 

1.32 

Heavy target 
fragments (Z>2) 

0 

.58 

.60 

.65 

.76 

.84 

Total 

7.08 

8.76 

8.28 

5.17 

2.71 

2.28 




Table III. SEU Contribution From Various Radiation 
Components Behind Several Shield Thicknesses for STS-51 
(upsets per day per computer) 


Aluminum thickness, 
g/cm 2 

0 

3 

5 

20 

55 

100 

Projectiles and 
projectile fragments 

1.30 

1.16 

1.08 

.62 

.17 

.04 

Light target 
fragments (Z= 1 ,2) 

0 

.36 

.39 

.51 

.72 

.89 

Heavy target 
fragments (Z>2) 

0 

.24 

.26 

.29 

.37 

.46 

Total 

1.30 

1.76 

1.73 

1.42 

1.26 

1.39 


Table IV. Comparison of Flight Data 
and Analyses for SEU Rates 
(upsets per day per computer) 



Flight Data 

Previous Analysis 

This work 

STS-51 

2.13 

1.33 

1.52 

STS-56 

6.05 

6.26 

5.85 


VI. CONCFUSIONS 

A more consistent set of predictions for SEU rates 
observed in two Shuttle flights for low and high inclination 
orbits is obtained in this study than prior analysis which does 
not contain target fragmentation. The current analysis uses 
HZETRN which has been compared favorably with flight data 
obtained by various types of detector responses not limited to 
microelectronic response. It is clear in the comparison of the 
HZETRN prediction with the D1 track data that in addition to 
a correct environmental model one must have an adequate 
nuclear data base and computational procedures to estimate 
the particle flux-energy spectra within the shield material. 
The uniqueness of the Langley transport code development 
lies in the requirement of laboratory validation with well 
defined ion beams and target materials [24,25,26]. This is 
especially important in the low and intermediate LET region 
where heavy target fragment contributions are modest 
(LET < 300 MeV/g/cm 2 ). At higher LET, the prediction of 
etchable tracks must account for limitations of the detector 
foils since many high LET tracks are lost in post exposure 
processing. Attention to such details are beginning to close 
the gap between flight data and prediction. A similar attention 
to detail has been required to understand the TEPC data as 
well. We likewise find that evaluation of target fragmentation 
effects in SEU must be approached by accounting for some of 
the limitations on the generation and collection of charge 
within the device. In the case of SEU, we not only treat the 
distribution of upset sensitivity within the device elements but 
must allow for the longitudinal and lateral characteristics of 


the track structure and Auger recombination within the track. 

Attention to such details will further close the gap between 

predictions and observations for the shuttle computer 

SEU data. 
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